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ABSTRACT 


The  equations  of  motion  are  derived  for  an  object  falling  under 
the  influence  of  gravity  and  aerodynamic  drag.  These  equations  are  used 
to  examine  the  response  of  the  ARCAS  parachute  to  hypothetical  wind 
profiles  in  the  region  extending  from  30  to  80  km. 

A  computational  scheme  for  determining  the  horizontal  wind  from 
the  observed  motion  of  any  wind  sensor  influenced  only  by  the  above  forces 
is  presented,  along  with  an  example  of  its  application  to  an  ARCAS  parachute 
flight. 
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I.  INTRODUCTION 


One  important  question  asked  by  atmospheric  scientists  involved  in 
the  analysis  of  rocketsonde  data  concerns  the  accuracy  of  the  wind  informa¬ 
tion.  At  present,  the  wind  information  available  to  the  analyst  (Data  Report 
of  the  Meteorological  Rocket  Network  Firings)  is  determined  directly  from 
the  motion  of  the  wind  sensor  without  any  corrections  applied.  In  addition, 
the  data  presentation  commonly  begins  at  some  height  considerably  below 
apogee  (the  highest  point  attained  by  the  sensor). 

The  purpose  of  this  report  is  twofold: 

1.  To  Investigate  the  response  of  a  wind  sensor  to  its 
environment. 

2.  To  demonstrate  the  feasibility  of  applying  corrections 
to  the  sensor  motion  in  order  to  obtain  the  real  wind, 
and  of  extending  the  wind  data  to  higher  altitudes. 

2,  EQUATIONS  OF  MOTION 

We  shall  assume  from  the  outset  tha*.  the  wind  sensor  moves  only 
in  response  to  the  force  of  gravity  and  aerodynamic  drag.: 

From  Newton's  second  law  of  motion,  we  can  write 
dVp  _  - 

m  =  mg  +  D  (1) 

where  m  is  the  mass  of  the  sensor,  dV  /dt  the  instantaneous  sensor 

P 

acceleration,  g  the  acceleration  due  to  graviiy  (assumed  constant)  and 
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D  the  drag  force. 

Figure  1  is  a  combination  force  and  velocity  vector  diagram  and 

illustrates  the  relationship  among  the  terms  in  (1)  and  their  relationship 

to  sensor  velocity  Vp  and  wind  velocity  V.  Also  shown  are  the  scalar 
■+ 

components  of  (V  -  Y  ). 

The  magnitude  of  the  drag  force  is  assumed  proportional  to  the 
square  of  the  velocity  of  the  air  relative  to  the  velocity  of  the  sensor  and 
can  be  expressed 

|B|  =  p  |V  -Vp|2,  (2) 

where  p  is  air  density,  CD  the  drag  coefficient,  and  A  a  characteristic 
area  of  the  sensor. 

From  Fig.  1,  the  drag  force  expanded  into  component  form  yields 
□  =  i  |D|  cos  a+j  |D|  cos  g  +  k  |D|  cos  y 

and,  further,  as 

D  =  T  |D|ur  |V- Vp|'l+T|B|  vr  |v- V?rJ+ k|5|  wr  IV-Vpl'1  (3) 

where  u^,  v^,  and  w^  are  the  components  of  air  speed  relative  to  the  sensor. 
Substituting  (2)  into  (3)  and  recognizing  that  | V -  Vp |  =  [u  2  +  v  2  +  w2  ]  ^ 
we  obtain 

D  =  (i  u  +  jv  +  kw)[u2+v2+w2l^pC_A/2.  (4) 

r  r  r  1  r  r  rJrD 

Substituting  (4)  into  (1)  and  expressing  the  result  in  component  form,  we  get 


COMBINATION  FORCE  ANO  VELOCITY  VECTOR  DIAGRAM 
APPLIED  TO  A  WIND  SENSOR  OF  MASS  m  WITH  VELOCITY 
\Tp  IN  WIND  FIELD  WITH  VELXITY  V  AND  ACTING 
UNDER  INFLUENCE  OF  GRAVITY 'irig  AND  DRAG  FORCE  IT 

FIG.  I. 
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(10) 
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It  is  significant  to  observe  that,  for  the  condition  of  zero  vertical  wind, 
the  error  in  the  horizontal  wind  components,  u^  and  v^,  can  be  determinea 
from  the  sensor  velocity  alone,  without  regard  to  the  physical  character¬ 
istics  of  the  medium  and  sensor  that  are  involved  in  K.  However,  since 
there  is  no  unique  sensor  velocity  profile,  this  approach  in  determining  u^ 
and  v^  can  be  utilized  in  real  situations  only.  An  example  of  employing  this 
approach  is  reserved  for  a  later  section,. 

In  order  to  permit  a  more  general  understanding  of  how  the  sensor  will 
respond  to  its  environment,  we  shall  therefore  postulate  the  state  of  the 


environment. 
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3.  QUALITATIVE  RESPONSE  CHARACTERISTICS 

Since  we  have  already  assumed  no  vertical  wind,  our  immediate 
concern  is  with  (11)  and  (12).  Because  they  are  similar,  we  shall  consider 
only  (11).  Upon  substituting  u^=  u  -  u^  (11)  becomes 


du 


u2  V2 


r  „  r,  ,  r  i  ‘  du 
dz  =  UrK  [1+  2  ]  +  dz 


(14) 


In  its  present  form  (14)  is  not  integrable  without  applying  some  con¬ 
ditions.  The  first  of  these  is  that  over  a  sufficiently  small  layer  of  integration 

2  2 

K  is  constant  and  the  second  is  that  u  «  w  . 

r  P 

If,  in  addition,  we  examine  the  particular  case  in  which  the  vertical 
wind  shear  is  zero,  then 


du 

~  =  K%  (15) 


and  its  solution  is 


“r(2)  =  uj>0)exp[ -K(zo-  z)],  (16) 

where  Uf(zo)  is  the  initial  wind  error  at  height  zq  and  z  <  zq.  Equation  (16) 
demonstrates  that  the  wind  error  decreases  exponentially  with  the  distance 
the  sensor  falls  through  the  atmosphere. 

Maintaining  the  same  prescribed  conditions  and  letting  the  vertical 
wind  shear  vary  sinusoidally,  (14)  becomes 


i 


da  2nA  , 

-t— ■  =  Ka  +  , — r— -  sin  [-^  (z  -  z)] 
dz  r  *-  L  1  L  '  o  11 


where  2tt  A  / L  is  the  amplitude  and  L  its  wavelength, 
m 


A  solution  for  steady- state  conditions  is  given  by  Middleton  and  Spilhaus 


(1953,  p.  66). 


^1  +  K2L2/4tt2 


:sintl7<zo'z)  ‘ tan"1  llJ’  (18) 


If  we  substitute  (u  -  u^)  for  and  denote  the  phase  angle  by  i|r,  we 


u  =  u  +  - 

/l  +  K2L2/4n2 


where  u  =  A  cost— (z  -z)]< 
m  L  o  1 


sinll7(Zo"Z>‘ 


This  result  shows  a  wind  speed  indication  of  reduced  amplitude  which 
lags  behind  the  real  wind  speed  fluctuations  in  the  environment.  As  expected, 
increasing  values  of  either  K  or  L  increase  the  amplitude  of  the  wind  speed 
indication  and  reduce  the  phase  angle.  Fig.  2  is  a  plot  of  (19). 


4.  QUANTITATIVE  RESPONSE  CHARACTERISTICS 

In  order  to  investigate  quantitatively  the  response  of  a  wind  sensor  to 
its  environment,  it  is  necessary  to  deal  with  the  physical  characteristics  of 
the  environment  and  sensor.  The  remainder  of  this  report  will  be  concerned 
with  the  ARC  AS  parachute. 
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The  parachute  is  launched  aboard  an  ARCAS  rocket,  ejected  in  the 
vicinity  of  80  km  and  tracked  by  radar.  Figure  3  has  been  constructed  from 
data  contained  in  a  report  by  the  Atlantic  Research  Corporation  (1962).  It 
should  be  noted  that  because  the  ejection  time  is  fixed  at  128  seconds,  the 
separation  of  the  sensor  from  the  carrier  can  (and  does)  occur  before,  at, 
or  after  apogee  (point  C).  For  this  reason,  and  because  of  yaw  and  pitch 
of  the  rocket,  the  direction  of  ejection  is  uncertain.  After  ejection,  all 
three  components  (motor,  nose  cone,  and  parachute)  are  tracked  until  the 
parachute  becomes  the  largest  radar  target.  No  information  is  in  general 
made  available  as  regards  the  "blossoming"  of  the  parachute,.  Winds  are 
normally  calculated  beginning  at  the  point  where  the  wind  sensor  is  deemed 
to  be  "wind  sensitive"  (point  D,  Fig.  3).  This  point  is  ascertained  subjec¬ 
tively  from  the  radar  track  by  noting  a 'discontinuity"  in  the  plotted  path. 

The  aerodynamic  performance  of  the  falling  parachute  is  dealt  with 
elsewhere  (Beyers,  et  al,  1962;  Whitlock  and  Murrow,  1964).  It  is  sufficient 
to  note  that  the  effective  cross-sectional  area  varies  with  descent  (by  virtue 
of  the  shroud  lines  twisting  and  the  pendulum  motion  induced  by  the  instrument 
package)  and  the  parachute  may  slip  and  glide  and  perhaps  collapse  altogether. 
These  effects  are  more  predominant  at  the  higher  altitudes.,  Since  we  are 
interested  in  altitudes  from  30  to  80  km,  the  variation  of  K  with  height  must 
be  taken  into  account.  We  shall  first  consider  the  vertical  variation  of  the 


drag  coefficient,  CD> 


SCHEME  OF  IDEALIZED  ARCAS  PARACHUTE  FLIGHT 


FIG.  3. 
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Wagner  (1964)  has  fitted  an  analytical  curve  to  an  average  fall  speed 
profile  based  on  many  trials  and  produced  the  expression 

w  =  -  Kj  exp  [  <K2  z  -  K3)\  -  K4]  ,  (20) 

where  the  K's  are  empirical  constants.  The  data  from  which  (20)  was 
deduced,  however,  were  taken  from  the  30-65  km  layer  only..  The  effect 
of  this  equation  above  65  km  can  be  determined  in  the  first  instance  by 
using  it  to  generate  a  profile. 

Let  us  imagine  that  we  have  released  the  parachute  at  80  km  into  a 
calm  atmosphere.  Then  ur  and  v^  are  zero  so  that  (13)  can  be  solved  for 
Cp  resulting  in 

,  dw 

CD  =  (2m/pA)  [-  -^4  -L-]  (21) 

p  w 

P 

The  Cp  profile  has  be  >n  computed  after  substituting  (20)  into  (21)  for 
Wp,  using  the  U.  S.  Standard  Atmosphere  density  (Duberg,  et  al,  1962), 

A  =  164045  cm^i  and  m  =  2330  gm..  In  Fig.  4  (curve  A)  the  drag  coefficient 
continually  increases  above  60  km.  We  see  no  apparent  justification  for  this 
on  physical  grounds.  Rather,  the  increase  is  a  consequence  of  extrapolating 
Wagner's  w^  profile  to  heights  where  its  application  was  not  intended.  As 
(20)  indicates,  w^  increases  without  bound  as  z  increases  This,  of  course, 
is  unrealistic  because  the  fall  speed  must  be  zero  at  apogee. 


A  —  BASED  ON  WAGNER'S  ( 1964)  FALL 
SPEED  PROFILE  AND  EXTENDED 
FROM  65  Km  TO  80  km 

B  —  PROPOSED 


PLOT  OF  HEIGHT  VS.  DRAG  COEFFICIENT 
FIG.  4. 


V 


Co 
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On  the  other  hand,  the  profile  should  be  acceptable  below  60  km 
because  it  is  based  on  empirical  w^  data.  In  order  to  obtain  a  plausible 
0^  profile  for  the  entire  30-80  km  region,  we  shall  assume  simply  that 
0^  varies  linearly  with  height  and  is  strongly  weighted  by  the  data  below 
60  km. 

The  expression  for  the  drag  coefficient  is 


C  =  1  -  0.  00533z  (22) 

where  z  is  in  kilometers.  Eq,'  (22)  is  also  shown  in  Fig.-  4  (curve  B). 

We  now  have  sufficient  information  to  derive  a  new  w^  profile  that 
will  be  more  consistent  with  expectations  above  60  km.  To  accomplish 
this,  we  need  first  to  integrate  (11)  and  (13).  The  integration  will  be  per¬ 
formed  over  a  layer  500  m  thick  in  which  we  shall  assume  a  mean  value  of 

K  and  the  two  terms  u  and  w  inside  the  radical  (let  v  =  0).  The  results 
r  p  r 

of  integrating  (11)  and  (13)  are,  respectively, 

u^  =  [  ur(zQ)  +  B  *  du/dz  ]  exp  [  -Bfz^  -  z. )]  -  S’*  du/dz  (23) 

and 

wr  =  {[  wr2(zo)-  g  B'1  j  exp[  -  2B(zo-  z)]  +  gB_1  j^2,  (24) 


_  _  2  \ 

B  =  K  [1  +  (ur/w  )  ] 


Equations  (23)  and  (24)  have  been  applied  to  successive  layers  500  m 
thick  commencing  at  80  km  and  proceeding  downward  to  30  km.  The 
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numerical  computational  procedure  for  each  layer  was  begun  by  computing 
B  from  an  initial  guess  of  u^  and  wf  .  was  computed  from  (22)  and  p 
from  the  Standard  Atmospheric  densuy.  Given  the  boundary  conditions 
ur(zQ)  and  w?  (zq),  and  the  mean  shear  du/dz,  u^  and  wf  were  computed 
for  the  bottom  of  the  layer  in  question.  The  values  of  Uf  and  wf  were  used 
to  obtain  new  values  of  u^  and  wy  and,  then,  an  "improved"  estimate  of  B. 
This  value  of  B  was  substituted  into  (23)  and  (24)  to  arrive  at  improved 
estimates  of  u^  and  w^.  The  iterative  procedure  continued  until  no  sig¬ 
nificant  change  in  B  resulted.  The  final  values  of  u^  and  were  used  as 
the  upper  boundary  conditions  for  the  next  layer  down. 

The  computed  u^  and  wy  profiles  for  zero  real  wind  ( u  (z )  =  0)  and 
the  boundary  conditions  u^  =  300  m  sec  *  and  =  Om  sec  *  at  80  km  are 
displayed  in  Fig.  5  along  with  Wagner's  fall  speed  profile  [Eq.  (20)].  A 
sensor  speed  of  300m  sec  *  imparted  by  the  rocket  at  apogee  is  not  unreal¬ 
istic.  The  similarity  of  the  two  fall  speed  profiles  below  60  km  is  a  direct 
consequence  of  utilizing  nearly  the  same  drag  coefficient  in  (24)  as  that  pre¬ 
dicted  from  Wagner's  fall  speed  equation  and  there  being  practically  zero 
wind  error.  Their  diversity  above  60  km  follows  from  opposite  reasons. 
Additionally,  Fig.  5  shows  that  the  parachute  must  fall  20  km  to  a  height 
of  60  km  before  the  wind  error  is  negligible. 

As  seen  in  (24)  and  (25)  the  fall  speed  is  dependent  on  the  wind  error 
u^.  The  significance  of  this  dependence  can  be  examined  by  postulating 
another  real  wind  profile  with  different  boundary  conditions  on  at  apogee 
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and  comparing  the  computed  wy  or  w^  profile  with  that  in  Fig.  5  (curve  B). 

The  assumed  real  wind  profile  is  shown  in  Fig.  6  (curve  A). 

The  ur  and  wf  profiles  were  computed  in  the  manner  previously  pre¬ 
scribed.  The  differences  between  the  u^  profiles  shown  in  Fig.  6  and 
between  the  w^  (or  w^)  profiles  shown  in  Fig.  7  are  the  result  of  making 
u^  -  0  at  apogee  in  one  case  and  uf  =  300m  sec  ^  in  the  other.  Figure  6 
illustrates,  as  did  Figure  5,  that  the  effect  on  the  speed  u  the  sensor  of  the 
impulse  imparted  to  it  by  the  rocket  is  almost  negligible  at  60  km.  In  Figure 

7  we  observe  that  despite  the  300  m  sec  *  difference  in  horizontal  speed  at 

apogee,  there  is  only  a  small  difference  in  the  fall  speed  profiles.  The 

maximum  difference  of  about  7%  occurs  near  72  km.  A  comparison  of 

curves  B  in  Figs.  7  and  5  (u(80)  =  300msec  *)  indicates  they  are  nearly  the 
r 

same;  the  maximum  difference  is  about  1%.  The  conclusion  ensues  that  the 
2  2 

effect  of  the  ratio  u^  /wy  on  the  fall  speed  is  slight  and  can  be  neglected. 

We  now  consider  one  additional  hypothetical  case  in  which  we  have  an 

8  km  vertical  wave  in  the  u-profile.  The  aforementioned  iteration  procedure 
among  (23),  (24),  and  (25)  produces  the  results  shown  in  Fig.  8..  The  wind 
error  at  the  top  of  the  profile  was  assumed  to  be  zero.  The  parachute  faii 
speed  profile  has  not  been  shown  since  it  is  virtually  identical  with  that  in 
Fig.  7  (curve  A).  Figure  3  reveals  that  even  aftei  falling  to  a  height  of 

60  km  the  phase  shift  is  45°  and  the  observed  amplitude  is  only  60%  of  the 
wind  amplitude.  The  sensor  must  fall  an  additional  20  km  before  it  is  almost 
completely  responsive  to  the  wind.  Of  course,  as  the  wavelength  decreases, 
the  response  of  the  parachute  also  decreases. 
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Figures  9  and  10  display,  respectively,  the  response  and  phase  angle 
dependence  on  both  wavelength  and  height.  In  Fig.  9  we  have  plotted  the 
response,  i.  e. ,  the  ratio  of  the  observed  amplitude  to  the  wind  amplitude, 
as  a  function  of  wavelength  at  various  heights.  The  ratio  of  the  respective 

2  2  2  V 

amplitudes,  as  determined  from  Eq.  (19),  is  1  /[  1  +  (4rr  IK  L  )]  2,  At  any 
given  wavelength,  the  response  increases  with  decreasing  altitude.  At 
8  km  only  10%  of  the  amplitude  of  the  longest  wa\  ilength  (16  km)  is  detected 
whereas  at  30  km  more  than  99%  of  the  amplitude  of  all  wavelengths  greater 
than  1  km  is  detected  by  the  parachute. 

In  Fig.  10  we  see  that  at  any  given  height  the  phase  angle  decreases 
with  increasing  wavelength.  At  80  km  the  phase  shift  is  greater  than  80° 
for  all  wavelengths  shown  whereas  at  30  km  the  phase  shift  is  less  than  5° 
for  all  wavelengths. 

5.  APPLICATION  TO  REAL  WIND  DATA 

As  pointed  out  in  Section  2,  if  we  assume  no  vertical  wind,  we  can 
deduce  the  wind  structuie  from  the  observed  motion  of  the  falling  parachute;: 
that  is,  by  substituting  (7)  into  (5)  and  (6),  we  obtain 

du  I  dw 

u  =  v  (WP  "df" )/[“df  +  gl>  (26) 

and 

dv  /  dw 

v  =  V<WP  +gJ-  <27> 

We  shali  apply  (26)  and  (27)  in  finite-difference  form  to  an  actual 


FREQUENCY  RESPONSE  OF  ARCAS  PARACHUTE  vs 
WAVELENGTH  AT  VARIOUS  HEIGHTS 


FIG.  9 


21 


0  30  40  SO  '  60 

\(/  (DEG) 

HEIGHT  vs  PHASE  ANGLE  \jf  =  tan1 
VARIOUS  WAVELENGHTS 


FIG.  10 


23 


sounding  from  the  ARCAS  parachute.  The  data  (parachute  position,  velocity, 
and  acceleration)  are  contained  on  magnetic  tape  and  have  had  a  101-point 
filter  applied  to  the  various  parameters  (Webster,  1964).  The  frequency 
response  of  this  filter  can  be  deduced  from  the  numerical  operator  employed; 
however,  let  it  suffice  here  to  say  that  the  smoothing  effected  is  very  little, 
indeed,  even  at  the  high  frequency  end  of  the  spectrum..  Equation  (26)  can 
be  written  in  finite  difference  form  as 

u  (i+1)  -  u  (i-1)  /  w  (i+l)-w  (i-1) 

u(i)  =  up(i)  -  [wp(i)  ”E— ^ - ]/[-£-lSi -  +g]'  (28) 

where  i  denotes  the  current  value,  i-1  the  previous  value,  and  i+1  the  subse¬ 
quent  value.,  At  =  0.1  second  and  g  =  32,174  ft/sec2.  An  analogous  relation 
holds  for  the  v-component  (Eq.  (27) ). 

Because  the  sensor  falls  at  a  variable  speed  (about  4  times  as  fast  at 
60  km  as  at  40  km)  and  the  data  are  spaced  in  equal  time  increments,  a  sym¬ 
metrical  filter  applied  to  the  data  is  automatically  a  function  of  height. 

A  second  101-point  symmetrical  filter  has  been  designed  to  eliminate 
certain  undesired  waves  in  the  data.  The  theoretical  (continuous (frequency 
response  of  this  filter  is  given  by: 

R'(f)  =  exp(-l6f)2  (29) 

where  f  is  in  cycles  per  second.  The  corresponding  weight  function  provid¬ 
ing  this  response  is 

w(t)  =  (tt^2  /  16)e-ip(-TTt  /  16 )2, 


(30) 
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when  t  is  in  seconds. 

The  response  to  a  discrete  101-point  symmetrical  filter  is  given  by 

50 

R(f)rW(0)  +  2  £  W(j)cos  (2rgf),  (31) 

j=l 


where  the  W(j)  are  determined  from  (30)  subject  to  the  constraint* 


50 


’W(O)  +  2  £  W(j)  =  1. 
j=l 


Equations  (29)  and  (31)  are  plotted  in  Fig,.  11.  Fig.  12  displays  the  response 

as  a  function  of  wavelength  at  various  heights  assuming  the  fall  speed  profile 

in  Fig..  7  (curve  B).  We  note  that  for  waves  in  u  and  v  lesB  than  125  m  in 

P  P 

wavelength, less  than  10%  of  their  amplitude  is  retained  at  30  km,  whereas 
at  70  km  the  same  percentage  applies  to  wavelengths  less  than  2250m.. 

The  results  of  applying  (28)  and  its  counterpart  to  the  smoothed  u^, 
v  ,  and  Wp  data  for  the  flight  under  consideration  are  shown  in  Fig.  13., 
Below  about  58  km  the  difference  between  the  associated  horizontal  com¬ 
ponents  of  the  parachute  and  wind  speed  are  small.  Above  this  level  the 
differences  increase  to  a  height  of  about  72  km.  From  this  point  to  apogee 
(75  km)  the  computed  wind  speed  components  are  very  erratic.  Not  all  the 
computed  wind  speed  components  are  plotted  for  this  region  because  their 
magnitude  was  off  the  scale  of  the  graph. 


For  further  information  on  filtering  techniques  refer  to  Holloway  (1958). 


RESPONSE  AS  A  FUNCTION  OF  WAVELENGTH  FOR 
VARIOUS  HEIGHTS 


FIG.  12 
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An  examination  of  the  vertical  accelerations  of  the  parachute  in  the 

latter  region  reveal  that  they  fluctuate  rapidly  in  time  from  values  loss  than 

-2  -2 

g  to  values  greater  than  g;  in  one  instance  from  15  ft  sec  to  45  ft  sec  in 
one  second.  As  seen  in  Eqs.  (26)  and  (27),  when  the  vertical  acceleration  of 
the  parachute  approaches  g  the  correction  term  approaches  infinity..  Whether 
these  accelerations  are  real  or  are  a  consequence  of  the  inability  of  the 
radar  to  correctly  track  the  parachute  is  not  known.  If  the  accelerations  are 
real  then  it  may  be  that  there  are  other  forces  present  that  have  not  been 
accounted  for  in  Eq.  (1). 

The  response  features  exhibited  in  Fig.  13  are  consistent  with  those 
described  in  the  previous  two  sections  of  this  report,  at  least  up  to  about 
68  km.  It  is  quite  apparent  that  large  corrections  must  be  applied  to  the 
measured  winds  above  60  km  in  order  to  obtain  the  true  winds.  In  the 
example  considered,  the  corrections  varied  from  about  40  m  sec  to  about 
100  m  sec  ^  in  the  60-70  km  region. 

6.  SUMMARY 

This  report  has  dealt  with  the  response  of  a  wind  sensor  acting  only 
under  the  influence  of  aerodynamic  drag  and  gravity.  The  following  assump¬ 
tions  were  made:  zero  vertical  wind,  the  physical  characteristics  of  the 
sensor  and  environment  constant  with  height,  and  the  ratio  of  the  wind 
error  (difference  between  the  horizontal  wind  and  parachute  speeds)  to  the 
fall  speed  small  compared  to  unity.  Under  these  conditions,  it  was  found 
that  the  equation  for  the  response  of  the  wind  sensor  was  of  the  same  form 
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as  that  for  the  response  of  an  ordinary  mercury  thermometer.. 

Next,  the  equations  of  motion  (assuming  zero  vertical  wind)  were 
applied  to  the  ARCAS  parachute.  The  Standard  Atmosphere  density  end  a 
semi-empirical  drag  coefficient  profile  were  err  ployed  for  the  30  to  80  km 
region.  By  postulating  various  wind  profiles  for  this  layer  and  different 
ejection  speeds  at  apogee  (80  km)  the  following  results  were  established: 

1.  in  cu3es  of  zero  to  moderate  horizontal  wind  shear,  che 

parachute  must  fall  from  80  km  to  approximately  60  km 
before  it  becomes  essentially  "wind  sensitive"  (wind  error 
on  the  order  of  a  few  meters  per  second); 

2.  the  parachute  fall  speed  profiles  are  all  approximately  the  same; 

3.  the  effect  of  the  ratio  of  the  wind  error  to  the  fall  speed  on  the 

fall  speed  is  small. 

Last,  and  most  significant,  a  simple  numerical  procedure  was 
developed  from  the  equations  of  motion  (assuming  zero  vertical  wind)  that 
can  be  used  to  compute  the  wind  structure  from  the  observed  sensor  motion. 
This  procedure  involves  neither  the  physical  characteristics  of  the  sensor 
nor  those  of  the  environment.  Of  additional  significance  is  the  fact  that 
this  procedure  does  not  depend  on  the  power  of  the  velocity  term  in  the 
equation  for  the  drag  force. 

This  technique  was  applied  to  an  actual  ARCAS  flight  and  the  computed 
component  wind  speed  corrections  were  as  large  as  100 in  sec  ^  in  the  60-70 


km  region. 
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